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ABSTRACT 

The spectra of disc accreting neutron stars generally show complex curvature, 
and individual components from the disc, boundary layer and neutron star surface 
cannot be uniquely identified. Here we show that much of the confusion over the 
spectral form derives from inadequate approximations for comptonization and for the 
iron line. There is an intrinsic low energy cutoff in comptonised spectra at the seed 
photon energy. It is very important to model this correctly in neutron star systems as 
these have expected seed photon temperatures (from either the neutron star surface, 
inner disc or self-absorbed cyclotron) of ~ 1 keV, clearly within the observed X-ray 
energy band. There is also reflected continuum emission which must accompany the 
observed iron line, which distorts the higher energy spectrum. We illustrate these 
points by a reanalysis of the GINGA spectra of Cyg X-2 at all points along its Z 
track, and show that the spectrum can be well fit by models in which the low energy 
spectrum is dominated by the disc, while the higher energy spectrum is dominated by 
comptonised emission from the boundary layer, together with its reflected spectrum 
from a relativistically smeared, ionised disc. 

Key words: accretion, accretion discs, radiation mechanisms:thermal, binaries: close, 
stars: individual: Cyg X-2, stars: neutron, X-rays: stars 



1 INTRODUCTION 

Disc accreting neutron star systems are amongst some of 
the brightest X-ray sources in the sky. They can be split 
into two broad types namely atolls and Z sources (Hasinger 
& Van der Klis 1989). Observationally these differ in that 
they show distinct types of spectral variability and timing 
behaviour, which link more fundamentally to a difference 
in the accretion rate and probably also in magnetic field 
(Hasinger & Van der Klis 1989). The atolls have an accre- 
tion rate which is typically below ~ 10% of the Eddington 
limit, and have weak magnetic fields while the Z sources ac- 
crete at 50— 150 % of Eddington, and probably have stronger 
magnetic fields (e.g. the review by Van der Klis 1995). Both 
types generally show 2-20 keV spectra which exhibit pro- 
nounced curvature, probably due to Compton scattering of 
seed photons by electrons with temperature kT^ ~ 5 keV 
(e.g. White, Stella & Parmar 1988), except at low luminos- 
ity states which are only seen in the atoll systems, where the 
spectrum makes a transition to a power law (e.g. Mitsuda 
et al., 1989; van Paradijs & van der Klis 1994; Barret et al., 
2000). 

Here we concentrate on the X-ray spectra of the Z 



sources. A full description of the accretion fiow for these 
sources should take into account the interaction of the disc 
with the stellar surface, including the effects of radiation 
pressure and magnetic field in a relativistic potential. This 
is likely to be extremely complex, as shown by initial hy- 
drodynamic (Popham & Sunyaev 2001; Inogamov & Sun- 
yaev 1999) and magnetohydrodynamic simulations (Miller 
& Stone 1997). The boundary layer/inner disc emission can 
then illuminate the rest of the disc giving rise to a Comp- 
ton reflected spectrum and associated iron Ka line from the 
heated, ionised, upper layer of the disc (see e.g. the review 
by Fabian et al., 2000, with specific application to Z sources 
by Brandt & Matt 1994). Given this complexity it is perhaps 
unsurprising that the observed spectra from Z sources (and 
atolls) have proved difficult to fit (e.g. Mitsuda et al. 1984; 
White, Stella & Parmar 1988). This led to the widespread 
use of colour-colour diagrams as opposed to true spectral 
fitting. Here we show that previous difficulties were mainly 
caused by using poor approximations for the comptonised 
and reflected emission components. We include relativistic 
smearing (Fabian et al., 1989) of the disc emission as the 
source of the line broadening (which is otherwise very dif- 
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ficult to explain) and derive excellent fits to the GINGA 
spectra of Cyg X-2 at all points on its Z track. 



2 X-RAY EMISSION MODELS 

X-rays can arise from a variety of processes from accretion 
onto a low magnetic field neutron star. There can be com- 
ponents from the neutron star surface, the boundary layer 
and the X-ray illuminated disc, and these should be sub- 
ject to general and special relativistic effects. The thermal 
emission from the neutron star surface is probably buried be- 
neath a very optically thick boundary layer at the high mass 
accretion rates seen in the Z sources (Popham & Sunyaev 
2001), leaving only the disc and boundary layer as the major 
emission components. The energy release in the boundary 
layer heats the material, which then cools predominantly 
by comptonization of any seed photons. The seed photons 
can either be from the disc, or from the neutron star sur- 
face, or from self-absorbed cyclo/synchrotron photons from 
the magnetic field (Lamb 1989; Psahis, Lamb & Miller 1995 
but see Wardzinski & Zdziarski 1999). The resulting spec- 
trum is then either clearly comptonised, or quasi-blackbody 
if the optical depth is large enough to mostly thermalise 
the emission (saturated comptonization). Further Compton 
scattering, including bulk motion effects, can take place if 
the accreting material has a substantial scale height, so that 
the emission from the boundary layer has to travel through 
this infalling material before it can escape to the observer 
(Lamb 1989; Psaltis ct al., 1995). 

The boundary layer can directly illuminate some frac- 
tion of the disc. Part of the illuminating flux is reflected 
(Smale et al 1993; Brandt & Matt 1994), while the rest is 
thermalised, enhancing the disc emission (e.g. Czerny, Cz- 
erny & Grindlcy 1986). This heats the upper layers of the 
disc to the Compton temperature of around 1-2 keV. This 
exceeds the virial temperature at distances > 3 — 6 x 10^° cm 
for a IAMq neutron star, so the material escapes as a wind 
forming a corona around the outer accretion disc (Begel- 
man, McKee & Shields 1983). This accretion disc corona 
(ADC) is especially important in determining the observa- 
tional appearance of the highly inclined systems, where the 
intrinsic, central source is hidden by the disc but can be seen 
via scattering in this tenuous material above and below the 
disc (e.g. White, Nagase & Parmar 1995). The photo ionised 
ADC emits line and continuum radiation, and can also scat- 
ter some of the boundary layer/inner disc emission back onto 
the disc, increasing the illumination (e.g. Vrtilck, Sokcr & 
Raymond 1993). Even at somewhat smaller radii the X-ray 
illumination heats the upper layers of the discs, so it has an 
appreciably larger scale height than an unilluminated disc, 
although the stronger gravity means that it cannot escape to 
form a wind. The emission from this photo-ionised, heated 
layer consists of recombination features (line and continuum 
radiation) together with the reflection of some fraction of 
the incident X-rays from the disc (KaJhnan & White 1989; 
Vrtilek, Soker & Raymond 1993; Ko & Kallman 1994). 



3 SPECTRAL MODEL 

From the discussion above, the spectrum should consist of 
components from the disc and boundary layer, together with 
reflection from the disc and scattering/emission/absorption 
from the ADC. Generally progress is made by estimating 
which components dominate so as to form a simplified model 
which nonetheless describes most of the physical processes. 
This is where the problems arise as there is no consensus on 
which effects dominate. 

One approach to modelling the spectra (generally called 
the 'Eastern model' in the literature) assumes that an op- 
tically thick disc dominates at low energies, while a Comp- 
tonised boundary layer dominates at higher energies. This 
was developed by Mitsuda et al., (1984; 1989), and fits the 
data well with the addition of a broad iron Ka line (White, 
Peacock & Taylor 1985; Hirano et al., 1987; Asai et al., 
2000). 

However, spectral fitting is inherently non-unique, and 
another approach is to assume that the majority of the emis- 
sion at both high and low energies comes from a comptonised 
disc, while the boundary layer emits as a quasi-blackbody 
in the middle of this energy range (White ct al., 1988). Such 
models (generally termed 'Western' models) can fit the data 
equally well, but again require the inclusion of a broad iron 
line component (White et al., 1988; Mitsuda et al., 1989; 
Hasinger et al., 1990). 

Here, we follow the Eastern approach, as it seems most 
likely that the boundary layer is strongly Comptonised 
(Popham & Sunyaev 2001). We assume that the low energy 
X-ray continuum is dominated by the disc, but the disc spec- 
trum is not well known especially at accretion rates close to 
Eddington where radiation pressure must be important. The 
radial temperature distribution depends on whether the vis- 
cosity is proportional to the gas pressure or the total (gas 
plus radiation) pressure (Stella & Rosner 1984). Electron 
scattering in the atmosphere is also expected to be impor- 
tant, leading to weak Comptonization of the disc emission 
(Shakura & Sunyaev 1973), and the intrinsic disc spectrum 
should be further broadened by relativistic effects (Ebisawa 
et al., 1991). Detailed calculations of these distortions show 
that with a restricted bandpass, where only the high energy 
tail of the disc spectrum is seen, then the shape of the spec- 
trum can be well approximated by a simple multi-colour disc 
model with a colour temperature correction (Ebisawa et al., 
1991; Shimura & Takahara 1995; Merloni ct al. 2000). Thus 
we approximate the disc spectrum by the diskbb model in 
XSPEC, but also test that this assumption does not strongly 
bias our results by using other disc spectral models. 

We assume that the seed photons for the comptoniza- 
tion in the boundary layer come from the neutron star sur- 
face as well as the inner disc, so they can be at higher ener- 
gies than the observed disc emission. The typical tempera- 
tures expected from close to the neutron star at high mass 
accretion rates are ^ 1 keV, i.e. within the observed X-ray 
bandpass. Simple analytic approximations (e.g. the compst 
model or a cutoff power law) for the Compton scattered 
spectrum are only valid where the seed photon energy is 
very much lower than that of the hot electrons (e.g. Hua & 
Titarchuk 1995). Explicit seed photon comptonization mod- 
els are required (such as the comptt model in XSPEC), and 
these give spectra with significantly different curvature close 
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to the seed photon energies. Such models are starting to be 
used to fit broad band (generally SAX) data from Z and 
atoll systems (Guainazzi et al., 1998; Piraino et al 1999; 
in't Zand et al., 1999; di Salvo et al., 2000; 2001). However, 
the assumption in comptt is that the seed photons have a 
Wien distribution. Here instead we assume that the seed 
photons have a blackbody distribition and use an analytic 
approximation to the Comptonisation based on the solution 
of Kompaneets equation (Zdziarski, Johnson & Magdziarz 
1996), as described in Zycki, Done & Smith (1999). 

This comptonising cloud around the neutron star can 
illuminate the accretion disc, ionising it and producing a 
characteristic reflected component and iron Kq line. This 
was first detected in Active Galactic Nuclei, and extensively 
studied both in supermassive and in in Galactic Black Holes 
(see e.g. Fabian et al., 2000 and references therein). However, 
its application to spectra from neutron star systems is much 
less widespread. The atoll systems at low mass accretion 
rate have power law spectra in the 2-20 keV range, and 
reflection is sometimes included to fit these data (Yoshida et 
al., 1993; Piraino et al 1999; Barret et al., 2000). However, 
the curving spectra from higher accretion rate atolls and 
Z sources have not been fit with reflection models despite 
observational (Smale et al., 1993) and theoretical (Brandt 
& Matt 1994) motivation. 

Another issue of importance in modelling the reflected 
spectrum is relativistic smearing. The disc extends down 
close to the surface of the neutron star, so reflection from 
the inner disc with relativistic distortions on both the line 
and reflected continuum should be important (Fabian et al 
1989). This has been seen in AGN (Tanaka et al., 1995) and 
Galactic Black Hole Candidates (Zycki, Done & Smith 1997; 
1998; 1999; Done & Zycki 1999) but again has not been 
studied in high mass accretion rate neutron star systems. 
Instead the observed strong broad line at 6-7 keV (White et 
al., 1985; Hirano et al., 1987; Asai et al. 2000) is generally 
interpreted as arising from the ADC and illuminated outer 
disc (Kallman & White 1989; Raymond 1993), ahhough it 
is very difficult to make both the observed line intensity and 
width from such a corona (see e.g. the discussion in Smale 
et al., 1993). 

We calculate the reflected emission from the thcomp 
continuum as described in Zycki et al., (1999). This cal- 
culates the self consistent line emission as well as the angle 
dependent reflected continuum for any ionization state and 
iron abundance (all the other elements are assumed to have 
solar abundance: Morrison & McCammon 1983). Relativis- 
tic smearing (including light bending) is then applied to this 
total (line plus continuum) spectrum (Fabian et al., 1989). 
The resultant spectrum depends on the inclination, which 
we fix at 60° (Orosz & Kuulkers 1999). 

We assume that emission/absorption/scattering from 
the ADC above the outer disc is negligible. This is probably 
not the case for the low energy spectrum, as line emission 
is seen at ~ 1 keV in several LMXB's (e.g. Bautista et al., 
1998 and references therein) including Cyg X-2 (Kuulkers 
et al., 1997). However, here we are concerned with the 2-20 
keV spectrum, and the iron K line associated with the ADC 
is probably small as the K/L line ratio is low for the steep 
spectra characteristic of the Z sources (Kallman 1995). 




Soft Hardness 



Figure 1. A colour - colour plot showing the Z track of Cyg X— 2 
in the June 1988 data. The overlaid boxes show the data used to 
make the individual spectra 



Table 1. GINGA observations of Cyg X-2 in June 1988 



UT Start 



UT End 



88-Jun-ll 05:25:54 
88-Jun-12 22:30:26 
88-Jun-14 04:00:34 



88-Jun-lO 12:35:00 
88-Jun-ll 13:30:00 
88-Jun-13 13:15:00 

4 OBSERVATIONS 
4.1 GINGA data 

We use the Cyg X-2 Ginga data, re-extracting them from 
the original FRF's using the Ginga data reduction software 
at Leicester. Observations were performed in MPC2, MPC3 
and PC modes (see Hasinger et al., 1990); however, in this 
paper, we use only the MPC2 mode data in order to per- 
form a detailed spectral analysis. Data were rejected from 
periods when the SUD rate {i.e. the LAC count rate above 
24 keV) exceeded 8 counts/sec/detector, the magnetic cut- 
off rigidity was outside the range 7-20 GeV/c, the angle be- 
tween the LAC pointing direction and the Earth's horizon 
was outside the range 6-120 degrees, and the SOL2 rate (of 
the solid-state electron detector on board Ginga). Unfortu- 
nately, the intensity of Cyg X-2 during the observations was 
such that the register storing the detected count rate over- 
flowed occasionally before telemetry. This affected only data 
accumulated at medium bit-rate, and could be corrected af- 
ter examining each PHA spectra. The data were corrected 
for deadtime using the following formulae: 



2 



1 — jJ-DET X r 



(1 - A*DET X r)2 



where, A and (t^ are the true count rate and variance, fi 
and CT^ are the observed count rate and variance, /^det is 
the observed count rate for one detector integrated over all 
PHA channels, and r is the deadtime {— 1.8 x 10~* seconds). 
Background subtraction was performed using off-source data 
taken in a contemporary 3-4 month period, and correcting 
for the 37-day periodicity of the Ginga orbit {i.e. the "uni- 
versal" method of Hayashida et al., 1989; Williams et al., 
1992). The spectrum of Cyg X-2 is sufficiently soft as to 



© 0000 RAS, MNRAS 000, 000-000 



4 C. Done, P.T. Zycki and D. A. Smith 




channel energy (keV) channel energy (keV) 



Figure 2. The two panels show spectral fits and residuals to 02 (horizontal branch: left hand panel) and 08 (flaring branch: right hand 
panel) spectra using the models detailed in Table ^. The low energy component shown by the dashed line is the disc blackbody, the dotted 
line is the comptonised continuum spectrum, the dashed-dotted line is its reflection from the disc while the solid line shows the total 
spectrum. The softening of the high energy emission is clearly the major difference between the horizontal and flaring branch spectra. 



not contaminate the SUD count rate, unlike other X-ray bi- 
naries such as Cyg X-1. Finally, we note the the Galactic 
emission estimated from a couple of off-source regions near 
to Cyg X-2 is negligible. 

The background subtracted data were rebinned on a 
40 seconds timescale. The soft hardness (4.5-9.0 keV/1.1- 
4.5 keV) versus hard hardness (14.2-9.0 keV/4.5-9.0 keV) 
ratio for the June and October datasets are shown in Figure 
|l| These are identical to those presented by Hasinger et al., 
1990. In the Hasinger et al., analysis, spectra were extracted 
at various intervals during the observations; to avoid spectral 
variations, these spectra were of short exposure. Here we 
extract spectra from each box in Figure |l|, thus enabling us 
to obtain much longer exposures hut avoiding any significant 
spectral variations and so a much more detailed analysis can 
be made. For the remainder of this paper, we refer to the 
spectra as 01-09, respectively. 



5 SPECTRAL ANALYSIS 

The interstellar column density to Cyg X-2 determined 
through visual extinction to the secondary star and the rela- 
tive size of the X-ray dust scattering halo is 2.2 x 10^^ cm~'^ 
(Predehl & Schmitt 1995). This is consistent with the col- 
umn seen from ROSAT, BBXRT, ASCA and SAX spectral 
analysis so we fix this in all the following fits. 

Results from fitting models based on the ideas devel- 



oped in Section 3 are given in Table g. We obtain excellent 
fits to all the individual spectra (total = 136.8/180) with 
optically thick disc emission and a Comptonised boundary 
layer together with its reflection from an ionized disc. The 
spectral fits for two different positions on the Z track (02 
and 08) are shown in Figure ^ 

Previous papers have often modelled the iron spectral 
features using a broad line. Table |^ shows the results ob- 
tained by replacing the more physically motivated refiected 
spectrum with an ad hoc broad iron line. The line energy and 
intrinsic width are constrained to be between 6.4-7 keV, and 
less than a < 0.42 (corresponding to a FWHM< 1 keV), re- 
spectively, for stability. The fit is adequate (although the 
model fails to describe the spectra on the normal branch - 
fiaring branch transition) with xt = 211.9/180. An F test 
shows that refiection is 99.98 per cent more probable to be 
the correct description of the spectral features than a broad 
Gaussian iron line. 

Figure ^ shows how the derived continuum parame- 
ters evolve along the Z track using the best fitting reflec- 
tion model for the spectral features (Table ^). All fluxes are 
derived by integrating the relevant unabsorbed model spec- 
trum over 0.01-100 keV. The most significant change in the 
spectrum is the factor 2 increase in disc flux between spec- 
tra 4 and 5. However, despite this increase in accretion rate 
through the disc, the flux in the boundary layer stays con- 
stant to within 10 per cent, although its spectrum softens 
(seen both by an increase in spectral index and decrease in 
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Figure 3. The derived continuum parameters for the spectral fits 
to the Eastern model with full boundary layer comptonization 
and reflection (diskbb+thcomp+ref lection) given in Table ^ as a 
function of position along the Z track. All temperatures are given 
in keV, and all fluxes are in units of 10~* ergs cm~^ s~^. The 
only statistically significant changes are that the disc blackbody 
flux increases while the comptonised spectrum softens. 



electron temperature). Thus the the ratio of flux dissipated 
in the boundary layer also makes an abrupt change between 
spectra 04 and 05, from ~ 3 to ~ 1.6. 

Figure ^ shows the evolution of the soft component in 
more detail. The disc temperature is clearly consistent with 
staying constant as the flux in the soft component increases 
(upper panel), which would imply that the optically thick 
boundary layer extends progressively over more and more 
of the inner disc, so the the observed inner disc radius in- 
creases. However, the uncertainties are large and cannot ex- 
clude the Tdiac oc fdisc dependence expected for increasing 
accretion rate through a constant radius flow. The lower 
panel of flgure |^ shows that the seed photons for the comp- 
ton scattered component have a temperature which is again 
consistent with a constant but is signiflcantly higher than 
the temperature of the observed soft component from the 
disc. This implies that the disc photons we observe are not 
the source of seed photons, and again, the most probable ge- 
ometry is one where the optically thick comptonising cloud 
surrounds most of the neutron star surface and also perhaps 
the innermost disc. 

The reflection parameters are shown in Figure ^. These 
are poorly determined since reflection is a curving spectral 
component and so is difficult to disentangle from the other 
two curving spectral components from the disc and comp- 



tonised boundary layer. The spectral features help, but the 
relative contrast of these can be dramatically changed by 
ionization effects (and by relativistic smearing). All the pa- 
rameters are consistent with a constant value {fl/2n = 0.16, 
^ = 330 and Rin = 120-Rg). However, the parameters are 
highly correlated, and together they are not consistent with 
these constant values. Fixing the reflected fraction (as ap- 
propriate if the geometry remains the same) at the mean 
value of Q/2n — 0.16 requires a signiflcant change in the 
ionizationparameter, as shown by the offset open symbols 
in Figure |5[ Physically, what is happening is that the spec- 
tral features increase in intensity along the Z track, so either 
the solid angle subtended by the reflector increases or its 
ionization state increases. Since the solid angle here is flxed, 
then the ionization state increases. Alternatively of course, 
it could be that the ionization state remains flxed and the 
geometry changes or that both change. Unfortunately, the 
broadening of the spectral features is still largely uncon- 
strained, these data cannot determine the inner disc radius 
from relativistic smearing. 



6 ALTERNATIVE SPECTRAL MODELS 

While this model gives a good description of the observed 
spectrum, the underlying continuum form is not unique. 
Even within the same physical picture of disc accretion and a 
comptonizing boundary layer there are different approaches 
to modelling the data. There are several approximate forms 
for the Comptonised spectrum. Our model, which is very 
similar to the comptt model in XSPEC (Titarchuk 1994), in- 
cludes the spectral curvature expected from having the seed 
photons within the observed bandpass. Previous papers have 
often used the compst model (Sunyaev & Titarchuck 1980) 
or cutoff power law for the comptonised flux, but this is ap- 
plicable only where the seed photons for the Comptonization 
are at much lower energies than the observed bandpass. This 
is plainly not the case for the Z sources (nor for any neutron 
star accreting at more than 0.1 per cent of Eddington). 

In Table ^ we flt the phenomenological model of a 
disc blackbody, comptonised continuum without seed pho- 
tons and broad iron line) to the data, and obtain very dif- 
ferent results (compare to Table |^ where the comptonised 
model includes seed photon curvature). Firstly, the flts are 
much worse using the compst approximation - typically by 
Ax^ ~ 10 larger in each spectrum for one less model pa- 
rameter (the seed photon energy). More importantly, the 
derived model parameters are very different. Figure com- 
pares the continuum luminosities derived from the compst 
model (Table to those derived from the thcomp model 
(Table m) where the curvature in the comptonised spectrum 
near the seed photon energy is included. The disc tempera- 
ture is much higher in the compst fits, as is its flux, but its 
implied radius is much smaller. 

The mean radius of the inner disc derived from the 
thcomp approximation (including seed photon curvature) is 
6OI4Q km, compared with 9.2 ± 0.2 km from compst (from 
the mean normalisation of the diskbb components, assum- 
ing that cos{inclination) / D^Q ~ 1 as appropriate for the 
best estimates of distance in units of 10 kpc, Dio ~ 0.72, 
and inclination of 60°: Orosz & Kuulkers 1999). Both are 
of order the neutron star radius, especially given the many 
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Tdisk (kev) 

Figure 4. The soft spectral components from the Eastern 
model with full boundary layer comptonization and reflection 
(diskbb+thcomp+ref lection). The upper panel shows disc tem- 
perature versus disc flux from the flts in table ^. The disc temper- 
ature is consistent with a constant as the flux vaxies by a factor 
2, which would imply that the disc inner radius was increasing. 
However, the uncertainties are large and the data cannot rule out 
the expected flux-temperature dependence for a constant radius 
flow, shown by the solid curve. The lower panel shows the seed 
photon temperature as a function of the disc temperature. The 
seed photons are plainly significantly hotter than the observed 
disc temperature (the solid line shows T^eed = Tdiac)} showing 
that the observed soft component is not the major source of seed 
photons for the comptonizing cloud. 



systematic effects inherent in the diskbb radius (neglect of 
colour temperature correction, relativistic effects and stress 
free inner boundary condition), but the apparent emitting 
area is suppressed by a factor of 7 by using an incorrect 
spectral form for the comptonised spectrum. Neglecting the 
turnover at the seed photon energy means that the contri- 
bution of the comptonised emission is overestimated at low 
energies, so the remaining flux to be associated with the disc 
emission is underestimated. Previous comparisons the West- 
ern and Eastern models have cited the very small emitting 
area of the disc as a problem for the Eastern models (e.g. 
Church & Balucinska-Church 2001), but here it is clear that 
the small area can be an artifact of incorrect spectral mod- 
elling of the comptonised component. 

By contrast, changing the description of the disc emis- 
sion does not dramatically change the inferred parameters 
for fluxes from the disc and boundary layer components. The 
diskbb disc blackbody model approximates the disc spec- 
trum neglecting the stress-free inner boundary condition. 



w o 



H h 
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Spectrum number (proportional to S^) 

Figure 5. The derived reflected parameters from the Eastern 
model with full boundary layer comptonization and reflection 
(diskbb+thcomp+ref lection) given in Table ^ as a function of 
position along the Z track. All the parameters are statistically 
consistent with constant values due to the large uncertainties, 
although this masks subtle spectral changes. Assuming that the 
solid angle subtended by the reflector is constant then there are 
significant changes in the ionization of the reflected spectrum. 
These are shown by the open symbols which are offset by 0.1 in 
spectral number. 



relativistic corrections and any energy input from illumina- 
tion, and assumes that the disc spectrum at any radius is 
a true blackbody (Mitsuda et al 1984). The diskm model in 
XSPEC (Stella & Rosner 1984) tackles two of these assump- 
tions, in that it uses a modified blackbody spectrum to de- 
scribe the disc spectrum where the opacity is dominated by 
electron scattering (Shakura & Sunyaev 1973), and includes 
the inner boundary condition. The spectrum is then some- 
what broader than the diskbb model, but gives results which 
are more or less the same as those for the diskbb model 
in the 2-20 keV range in terms of the total disc fiux, and 
concomitant comptonised continuum parameters. We fix the 
diskm spectrum normalisation — cos{inclination) / D^q at 
unity (Orosz & Kuulkers 1999), alpha viscosity at 0.1 and 
mass of the neutron star at 1.4M0. Using the diskm model 
with the thcomp continuum and broad iron line to model 
the spectrum gives an adequate fit (total xt = 224.1/180), 
but replacing the broad Gaussian line by the refiected spec- 
trum again gives a significant improvement in the fit (total 
xl = 139.5/180). 

The differences between diskm and diskbb are quite 
subtle. However, if there is substantial comptonization of 
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Spectrum Number (proportional to S^) 

Figure 6. Continuum parameters derived from the Eastern 
model assuming approximate boundary layer comptonization 
(compst, dark grey, Table ^) compared to those derived from mod- 
els which include the seed photon curvature (thcomp, light grey, 
Table ^). Both sets of fits have the disk spectrum modelled by 
diskbb and reflection is approximated by a gaussian line. The 
disc temperature is given in keV, and all fluxes are in units of 
10~* ergs cm~^ s~^. Plainly there is a dramatic difference in the 
inferred continuum depending on whether the comptonised spec- 
trum includes the effect of the seed photons. Given that the seed 
photon are expected to be within the observed bandpass then it 
is clear that approximations to the comptonised spectrum which 
do not include this can be severely misleading. 



o . 
«^ O 




Spectrum Number (proportional to S^) 

Figure 7. Continuum parameters derived from the Western 
model assuming approximate disc comptonization (compst, dark 
grey. Table ^) compared to those derived from models which in- 
clude the seed photon curvature (thcomp, light grey. Table 
Both sets of fits have the thermalized boundary layer modelled 
by bbody and reflection is approximated by a gaussian line. The 
approximate comptonisation model requires a blackbody at ~ 1 
keV as appropriate for the boundary layer, but its luminosity is 
far below that of the disc. There is a substantial drop in the de- 
rived 'boundary layer' blackbody temperature to ~ 0.3 keV when 
the comptonised disc is modelled by the more physically realistic 
comptonization which includes the curvature expected at energies 
close to the seed photons. The blackbody temperature is given in 
keV, and all fluxes are in units of 10~* r--^~'^ 



the disc emission then the disc spectrum can change dra- 
matically. Such strong comptonization is sometimes seen in 
the Galactic Black Holes at high mass accretion rates (Zycki, 
Done & Smith 2001). This is the basis of the 'Western' ap- 
proach, in which it is the disc emission which is strongly 
comptonised, while the boundary layer is a (more or less) 
simple blackbody (White, Stella & Parmar 1988, see Section 
2). Table ^ shows results from this with the commonly used 
compst model, together with the phenomenological broad 
iron line to model the spectral features. This is a good fit 
to the spectra (total = 171.5/189), and has a blackbody 
boundary layer temperature of ~ 1 keV, but the ratio of 
the boundary layer to disc flux is much lower than that pre- 
dicted for a non-rotating neutron star. This was noted by 
White et al., (1988) and suggested as evidence for equilib- 
rium (fast!) neutron star spin periods. However, this result 
is not robust. The compst model assumes that the seed pho- 
tons are far below the observed bandpass, so the spectrum 
is a power law below energies corresponding to the electron 
temperature. But the seed photons must be in the observed 



energy range, since the assumption here is that it is the disc 
emission itself which is being comptonised. Using the thcomp 
model for the comptonised disc, so that the spectral curva- 
ture expected from the seed photons is included, gives very 
different results (Table ^). Figure ^ shows a comparison of 
the derived blackbody boundary layer temperature and flux 
for these two different disc comptonization models. Plainly 
the continuum parameters vary substantially, and again, we 
stress that using compst outside of the range where its ap- 
proximations are valid can give misleading results. 

Figure ^ compares the derived spectra for datasets 02 
and 08 using the black body boundary layer plus comp- 
tonised disc approach. The upper panel shows the fits where 
the seed photon curvature is neglected, which is the standard 
approach in the literature, while the lower panel shows how 
the fits change when the seed photon curvature is included 
in the models. It is clear that the blackbody boundary layer 
temperature changes dramatically, from ~ 1 keV derived 
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Figure 8. Spectra derived using tlie Western model, assuming that the disc emission is strongly comptouiscd while the boundary layer 
is mostly thermalised. The upper panel shows the derived fits to 02 (horizontal branch) and 08 (flaring branch) spectra using the compst 
model for the comptonised disc. This does not have the seed photon curvature included, and the data then require a high temperature 
~ 1 keV blackbody component to fit the observed spectra. The lower panel shows the same data deconvolved using a model which does 
include the seed photon curvature in the comptonised spectrum. This dramatically changes the derived boundary layer temperature, 
reducing it to ~ 0.3 kcV. Comptonization approximations which do not take into account the seed photons can give highly misleading 
results for the boundary layer emission. 
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Table 2. GINGA data fit to the Eastern model with comptonization of the boundary layer including seed photon curvature and full 
reflection (diskbb + thcomp + reflection) 



data fcTdiskbb keV /diskbb 



kT, keV fcT,eed keV 



C/27r 



01 
02 
03 
04 
05 
06 
07 
08 
09 



n cro + 0.30 
U.OCi_Q 20 
,+0.31 
-0.13 
1+0.90 
-0.15 
S+0.99 
-0.18 
1+0.15 
-0.13 
3+0.18 
-0.14 
™+0.18 
-0.14 
-+0.19 
-0.20 
^+0.28 
-0.23 



o.5i: 

0.61^ 
0.58^ 
0.4lj 
0.39^ 
0.36j 
0.35^ 
0.34^ 



0.41 
0.43 
0.44 
0.44 
0.68 
0.74 
0.89 
0.82 
0.64 



2.30 
2.32 
2.49 
2.53 
2. 

2.90 
2.77 
2.69 
2.54 



+0.22 
0.11 
+0.25 
0.08 
+0.32 
0.11 
+0.25 
0.17 
+0.12 
-0.16 
+0.30 
0.16 
+0.25 
0.14 
+0.21 
0.15 
0.19 
0.11 



3 S6+°-l'5 
'^■■"-'-0.14 

'-'•^"-0.10 
I +0.50 
-0.20 
1+0.36 
-0.26 
9 QS+0-30 
^■y^-0.19 

7+0.33 
-0.14 
2 74+0.26 
^- '^-0.14 
s+0.19 
-0.12 



3.41] 

3.31"' 



2.77"* 



2.66"* 



0.82^ 



0.801°:- 

1+0.21 
-0.07 

'^■''"-0.09 

0.87l°;f3 
0.78 ±0.06 



0.80 
0.80 
0.81 
0.80 



+0.08 
0.05 
+0.07 
0.05 
+0.07 
0.06 
+0.09 
0.05 



1.14 

1.27 
1.30 
1.26 
1.20 
1.26 
1.26 
1.28 
1.13 



,0+0.58 
'^■^'-'-0.09 

07+0-5'' 
'^■^'-0.03 
j+0.57 
-0.04 
,+0.15 
-0.06 
,+0.24 
-0.06 
-+0.10 
-0.09 

'^•^"-0.05 
,+0.18 
-0.05 
,+0.11 
-0.06 



0.08: 
0.11^ 
0.21^ 
0.25^ 



O.li: 
0.11^ 



no; 

670^ 
570 
9300 
1000 
240 
2100 



-9000 
-110 
18000 
670 
+20000 

570 
+60000 

9260 
+8000 
-900 
1600 



140 
13000 
2000 
16000 



I2OOI1110 

1.80t^-*4 X 



38 
42 
49 
49 
40 



+00 

32 



-36 
-00 
-43 



-43 
n,+240 
-25 

1000 ■ 



90 
92 
1000 



-985 

CXD 

75 
00 
80 
+00 



8.4/20 
9.0/20 
12.0/20 
14.6/20 
24.8/20 
15.6/20 
22.6/20 
12.4/20 
17.4/20 



" all fluxes are the unabsorbed model, integrated from 0.01-100 keV, in units of 10 ergs s cm 

Table 3. GINGA data fit to the Eastern model with comptonization of the boundary layer including seed photon curvature with approx- 
imate reflection (diskbb + thcomp + gau) 



data fcTdiskbb keV /diskbb 



fcTe keV fcT^eed keV f^. 



thcomp 



keV^ 



keV 



EW (eV) xV^^ 



01 


0.46 


02 


0.46 


03 


0.53 


04 


0.52 


05 


0.45 


06 


0.42 


07 


0.41 


08 


0.39 


09 


0.37 



+0.11 

-0.12 
+0.14 
-0.09 
+0.13 
-0.11 
+0.15 
-0.10 
+0.09 
-0.07 
+0.12 
-0.10 
+0.11 
-0.08 
+0.13 
-0.12 
+0.22 
-0.13 



0.47 
0.51 
0.48 
0.51 
0.73 
0.74 
0.78 
0.67 
0.59 



2.30 ±0.05 
2.35j 
2.53j 
2.67j 
3.30^ 
3.24j 
3.15j 
2.i 

2.72^ 



-+0.O6 

-0.04 
ci+0.08 
-0.07 
7+0.12 
-0.10 
-1+0.20 
-0.14 
^+0.25 
-0.20 
-+0.22 
" -0.14 
-,+0.17 
='-0.12 
,+0.17 
-0.09 



3.42 
3.34 
3.50 
3.58 
3.77 
3.53 
3.32 
3.00 
2.83 



+0.11 
0.10 
+0.11 

0.07 
+0.17 
0.13 
+0.30 
0.19 
+0.62 
0.34 
+0.70 
0.43 
+0.50 
0.34 
+0.27 
0.18 
+0.22 
0.12 



s+0.06 
-0.05 
j+0.05 
-0.04 

3+0.07 
-0.05 
™+0.05 
-0.03 
7+0.05 
-0.04 
7+0.05 
-0.03 
-+0.05 
-0.04 
-+0.06 
-0.04 



0.78: 
0.83^ 



0.89: 
0.86^ 
0.87j 
0.87j 
0.86j 
0.85^ 



1.20 
1.31 
1.36 
1.36 
1.24 
1.29 
1.31 
1.30 
1.20 



6.7 ±0.3 

6.8 ± 0.2 

6 9+0-1 
"■^-0.3 

7.0-0.2 

6.8l°- 

6.7 ±0.2 

6.8 ±0.2 
6.8 ±0.2 

7.0-0.2 



0.42_o.4i 

0.42-0.42 
0.42-0.42 
0.42-0.27 
0.42-0.08 
0.42-0.10 
0.42-0.14 
0.42-0.42 
0.42-0.42 



50 ± 18 
52-22 

^o+i'' 

0*^-20 

97 ±20 

looli 

+0.20 

25 
+ 15 
25 
+ 16 
23 



86 
60 
50 



11.0/20 
12.5/20 
15.1/20 
20.2/20 
47.7/20 
29.5/20 
38.8/20 
15.8/20 
21.3/20 



" all fluxes are the unabsorbed model, integrated from 0.01-100 keV, in units of 10" 
' Energy constrained to be between 6.4—7 keV 
Width constrained to be between 0-0.42 keV 



ergs s cm 



from the inapplicable compst model to ~ 0.3 keV, in models 
with the seed photon compton curvature (thcomp). 



7 DISCUSSION AND CONCLUSIONS 

The spectrum of Cyg X-2 has complex curvature which 
makes spectral fitting highly non-unique. It is generally 
agreed that the X-ray spectrum must consist of at least 
two components, the disc and the boundary layer, and that 
at least one of these (though more probably both) is comp- 
tonised. In addition, there are spectral features correspond- 
ing to ionised iron K, which are often modelled as a broad 
Gaussian line. We have shown that the derived spectral pa- 
rameters are crucially dependent on how the comptonised 
spectrum is modelled. The often used compst or cutoff power 
law models are physically inconsistent since the seed photons 
for the compton scattering must be in the observed X-ray 
range, whether they come from the disc or the boundary 
layer. The lack of curvature in compst or cutoff power laws 
at energies below the electron temperature then forces a 
blackbody at ~ 1 keV into the fit. This is purely an artifact 
of incorrect modelling of the comptonised spectrum and has 



no physical relation to the boundary layer or neutron star 
surface. 

When the comptonised spectrum is modelled including 
this curvature (using thcomp or comptt) then the spectrum 
is adequately fit either by assuming that the boundary layer 
is comptonised, and the softer emission is from the accre- 
tion disc (diskbb+thcomp+gau xi = 211.9/180) or that the 
disc spectrum is comptonised, and the boundary layer is a 
blackbody (bb+thcomp+gau xl = 207.3/180). The GINGA 
data cannot distinguish between these models statistically. 
However, there are some theoretical constraints on the ra- 
tio of the disc and boundary layer fluxes. In the 'Eastern 
model', where the boundary layer is comptonised, we derive 
a boundary layer luminosity which is about twice that of the 
luminosity of the disc. This is as expected in the general rel- 
ativistic models of the gravitational potential if the neutron 
star spin is not close to the Keplarian breakup frequency 
(Sunyaev & Shakura 1986). By contrast, in the 'Western 
model' where the disc is strongly comptonised, the ratio of 
the boundary layer blackbody luminosity to that from the 
disc is ~ 0.2 — 0.25, implying that the neutron star is signif- 
icantly spun up. Detailed calculations indicate that to get a 
ratio of boundary layer to disc flux this low requires that the 
spin frequency of the neutron star is > 1 kHz (Sibgatullin & 
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Table 4. GINGA data fit to the Eastern model with approximate comptonization of the boundary layer and approximate reflection (diskbb 
+ compst + gaussian) 

data fcTdi^kbb keV /diskbb I'T, (keV) f^^^^^^, Eunc kev| a keV | EW (eV) xV^^ 

01 l-63lo;o5 2.93 ±0.07 14.5+^-0 1.06 6.7 ± 0.2 0.42_o.42 50^11 10.4/21 

02 1.68 ±0.04 0.94 2.86I0 05 ^^-^-I'^s 0-82 ^-^to.l 0.42_o.42 471^4 13.0/21 

03 1.70 ±0.04 1.14 2.87^0.}? ^"^-"^-I'l 0-68 ^-^-o.l 0.42_o.42 48^" 16.9/21 

04 1.68 ±0.04 1.24 2.82+JJ;^^ 18.3+;^q 0.56 ^-^to.l 0.42_o.28 53+ 21.9/21 

05 1.51 ±0.02 1.50 2.44+JJ;g5 200+^5 0.16 ^-^to.l 0.42_o.o6 105 ± 20 50.6/21 

06 1.58lo;o5 1-53 2.55lo^J 200l^o 0-1'' ^-^-o.i 0.42_o.io 95 ± 20 41.0/21 

07 1.61^o;o2 1-55 2.58+Q |!|g 200+^^ 0.14 6.6 ± 0.2 0.42_o.i4 73+^2 59.7/21 

08 1.66lo;o2 1-49 2.59 ±0.09 2001^^5 0-1'' ^-^-o.i 0.42_o.42 50^20 47.6/21 

09 1.72j;o;o3 1-35 2.63lo;Jg 200lj^5 0-13 6.7 ±0.3 o+0-i2 22t\l 60.8/21 



" all fluxes are the unabsorbed model, integrated from 0.01-100 keV, in units of 10~® ergs cm~^ 
^ Energy constrained to be between 6.4—7 keV 
Width constrained to be between 0-0.42 keV 



Table 5. GINGA data fit to the Western model with approximate comptonization of the disc and approximate reflection (bbody 
+ gau). 



compst 



data fcTbb (keV) /b 



kT, (keV) 



compst 



ke' 



keV 



I EW (eV) xV'^ 



01 


1.16 ±0.03 


0.26 


2.95 ±0.05 


12.0 ±0.3 


1.86 


6.8 ±0.2 


0.2 ±0.2 


35 ± 17 


7.5/21 


02 


1.18 ±0.03 


0.34 


2.87 ±0.04 


12.4 ±0.3 


1.92 


6 9+°-i 
"•"-0.3 


S+^'i 

^■■^-0.3 


oQ+20 


9.7/21 


03 


1.19 ±0.03 


0.41 


2.84_(, ,,4 


12.2 ± 0.3 


2.02 


7.0-0.3 


0.4-0.4 


34 ±20 


17.0/21 


04 


1.16 ±0.03 


0.44 


r, 70 + 0. 08 

^■'^-0.06 


12.1 ± 0.4 


2.06 


7.0-0.2 


0.4-0.2 




19.4/21 


05 


1.09 ±0.02 


0.51 


2.55 ±0.07 


11.1 ± 0.4 


2.45 


6.8 ±0.1 


0.4-0.1 


80l?° 


32.2/21 


06 


1.09 ±0.03 


0.55 


2.46 ±0.10 


11.9 ± 0.6 


2.12 


6.7 ±0.2 


0.4-0.1 


90 ±25 


21.2/21 


07 


1.10 ±0.03 


0.55 


2.44 ±0.08 


12.2 ± 0.6 


2.01 


D.M_Q 2 


0.4-0.2 


70 ±20 


30.9/21 


08 


1.10 ±0.03 


0.52 


2.38 ±0.06 


13.1 ± 0.6 


1.73 


6.8 ±0.2 


0+0.4 


38l?» 


11.9/21 


09 


1.09 ±0.03 


0.46 


2.34 ±0.05 


14.0 ±0.6 


1.44 


7.0-0.2 


0.4-0.4 




21.7/21 



" all fluxes are the unabsorbed model, integrated from 0.01-100 keV, in units of 10~* ergs s~^ cm"-^ 
^ Energy constrained to be between 6.4-7 keV 
Width constrained to be between 0-0.42 kcV 



Table 6. GINGA data fit to the Western model with full comptonization of the disc and approximate reflection (bbody + thcomp + gau). 

data fcTbb (keV) /bb T (keV) fcTe fcT,,,d fthco^p ^Unc ke\f] a keV | EW (eV) xV^^ 

01 0.36^0.05 0.24 2.29lo;o4 3.42lg;J° 0.78 ± 0.05 1.19 6.8 ± 0.2 0.42 - 0.42 50^20 10.5/20 

02 0.38 ±0.06 0.25 2.36 ± 0.05 3.35l°;Jg 0.83 ± 0.04 1.31 6.8 ± 0.2 0.42- 0.42 50^20 11.9/20 

03 0.41 ±0.06 0.26 2.52 ±0.07 3.49lg;i2 0.88 ± 0.04 1.37 ^-^to.l 0.42- 0.42 48^20 14.4/20 

+0.28 n «Q + n n/i 1 ■^7 7 n „ „ n /lo „ .„ ^1 +18 





fa 

J thcomp 


Ainc keVg 


0.78 ± 0.05 


1.19 


6.8 ±0.2 


0.83 ± 0.04 


1.31 


6.8 ±0.2 


0.88 ±0.04 
0.89 ±0.04 


1.37 
1.37 


6 9+''-^ 

D.M_Q 2 

7.0-0.2 


0.86 ±0.03 
0.87 ± 0.04 


1.25 
1.29 


6.8l»- 
6.7 ±0.2 


0.87 ± 0.04 
0.86 ± 0.04 


1.31 
1.30 


6.8l»;? 
6.8 ±0.2 


0.86 ± 0.04 


1.43 


7.0-0.2 



b 



04 0.41 ±0.06 0.27 2.67 ±0.10 3.57Io 2o 0.89 ± 0.04 1.37 7.0-0.2 0.42- 0.42 5II20 19.5/20 

05 0.37 ±0.05 0.36 3.31+q;j4 3.79+q;34 0.86 ± 0.03 1.25 6.8+q;2 0.42-o.io 95 ± 20 46.7 /20 

06 0.35 ±0.07 0.35 3.25lo;2o 3.54lJ5 j| 0.87 ± 0.04 1.29 6.7 ± 0.2 0.42_o.io 100 ± 25 29.1/20 

07 0.34 ±0.06 0.36 3.15lo;^3 3.32lg:*6 0.87 ± 0.04 1.31 6.81°;^ 0.42_o.i4 85I25 38.3/20 

08 0.32 ±0.08 0.33 2.891° 2.99lo^4 0.86 ± 0.04 1.30 6.8 ± 0.2 0.42- 0.42 64^3^ 15.7/20 

09 0.34l°;5^ 0.23 2.7it°o^{^i 2.861° ijg 0.86 ± 0.04 1.43 7.0-0.2 0.42-0.42 45 ± 21 21.2/20 

all fluxes are the unabsorbed model, integrated from 0.01-100 kcV, in units of 10~** ergs s^^ cm^^ 
Energy constrained to be between 6.4—7 keV 
Width constrained to be between 0-0.42 keV 
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Sunyaev 2000). This seems excessive: all known millisecond 
pulsars (the descendants of the LMXB neutron stars) have 
spin frequencies < 700 Hz (e.g. Kramer et al., 1998). 

Thus it seems likely that the 'Eeistern' model (Mitsuda 
et al., 1984; 1989) is correct and that the boundary layer 
dominates the hard spectrum, while the disc dominates at 
softer energies. A standard objection to this approach, that 
in general the derived disc parameters require an emitting 
area which is far too small to be a true disc (e.g. Church 
& Balucinska Church 2001), is again purely an artifact of 
incorrect modelling of the curvature due to the seed photons 
in the comptonised continuum. 

Finally, it is clear that the residual spectral features at 
iron are considerably better described by a strongly ionised 
reflected spectrum than by a broad Gaussian line. 
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